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Abstract 

We consider additive codes over GF(4) that are self-dual with respect 
to the Hermitian trace inner product. Such codes have a well-known 
interpretation as quantum codes and correspond to isotropic systems. It 
has also been shown that these codes can be represented as graphs, and 
that two codes are equivalent if and only if the corresponding graphs are 
equivalent with respect to local complementation and graph isomorphism. 
We use these facts to classify all codes of length up to 12, where previously 
only all codes of length up to 9 were known. We also classify all extremal 
Type II codes of length 14. Finally, we find that the smallest Type I 
and Type II codes with trivial automorphism group have length 9 and 12, 
respectively. 
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1 Introduction 

An additive code, C, over GF(4) of length n is an additive subgroup of GF(4) n . 
C contains 2 k codewords for some < k < 2n, and can be defined by a k x n 
generator matrix, with entries from GF(4), whose rows span C additively. C is 
called an (n, 2 fe ) code. We denote GF(4) = {0,1 , lu,lo 2 }, where lu 2 = uj + 1. 
Conjugation of x G GF(4) is defined by x = x 2 . The trace map, Tr : GF(4) — > 
GF(2), is defined by Tr(a;) = x + x. The Hermitian trace inner product of two 
vectors over GF(4) of length n, u = (u\, U2, ■ ■ ■ , u n ) and v = (i>i, Wa, . . . , v n ), is 
given by 

n n 

u*v = Tr(u ■ v) = Tr(nj^) = y^tt^ + u\vj) (mod 2). (1) 

i=l i=l 

Note that u * v is also the number (modulo 2) of places where u and v have 
different non-zero values. We define the dual of the code C with respect to 
the Hermitian trace inner product, C 1 - = {u G GF(4)™ | u * c = for all c G 
C}. C is self-orthogonal if C C C . It has been shown that self-orthogonal 
additive codes over GF(4) can be used to represent quantum error- correcting 
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codes |4|. If C = C^, then C is self-dual and must be an (n, 2") code. Self- 
dual additive codes over GF(4) correspond to zero-dimensional quantum codes, 
which represent single quantum states. If the code has high minimum distance, 
the corresponding quantum state is highly entangled. 

The Hamming weight of it, denoted wt(u), is the number of nonzero com- 
ponents of u. The Hamming distance between u and v is wt(u — v). The 
minimum distance of the code C is the minimal Hamming distance between any 
two distinct codewords of C. Since C is an additive code, the minimum distance 
is also given by the smallest nonzero weight of any codeword in C. A code with 
minimum distance d is called an (n,2 k ,d) code. The weight distribution of the 
code C is the sequence (Aq, Ai, . . . , A n ), where Ai is the number of codewords 
of weight i. The weight enumerator of C is the polynomial 

n 

W{x,y) =Y,A i x n - i v i (2) 

i=a 

We distinguish between two types of self-dual additive codes over GF(4). A 
code is of Type II if all codewords have even weight, otherwise it is of Type I. 
It can be shown that a Type II code must have even length. Bounds on the 
minimum distance of self-dual codes were given by Rains and Sloane |l9t The- 
orem 33]. Let di be the minimum distance of a Type I code of length n. Then 
dj is upper-bounded by 

if n = (mod 6) 

if n = 5 (mod 6) (3) 
otherwise. 

There is a similar bound on du, the minimum distance of a Type II code of 
length n, 




d H < 2 



A code that meets the appropriate bound is called extremal. It can be shown 
that extremal Type II codes must have a unique weight enumerator. Rains and 
Sloane a l so used a linear programming bound, and showed that extremal 
codes do not exist for all lengths. For instance, there is no self-dual (13, 2 13 , 6) 
code. If a code has highest possible minimum distance, but is not extremal, it 
is called optimal. An interesting open problem is whether there exists a Type II 
(24, 2 24 , 10) code. 

A linear code, C, over GF(4) which is self-dual with respect to the Hermitian 
inner product, i.e., u-v = for all u,d£C, is also a self-dual additive code with 
respect to the Hermitian trace inner product. However, most of the self-dual 
additive codes are not linear. Only Type II codes can be linear, since self-dual 
linear codes over GF(4) must contain codewords of even weight only. It follows 
that the set of Hermitian self-dual linear codes over GF(4) is a subset of the set 
of Type II self-dual additive codes over GF(4). 

Example 1. The unique extremal (6, 2 6 ,4) code, also known as the Hexacode, 
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has a generator matrix 



/l 
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1 LO LO 
LO LO 2 LO 2 



This code has weight enumerator W(x, y) 



45 x y + 18y b . It is therefore 



of Type II, and it can be verified that it is also a linear code. 

Two self-dual additive codes over GF(4), C and C', are equivalent if and only 
if the codewords of C can be mapped onto the codewords of C by a map that 
preserves self-duality. Such a map must consist of a permutation of coordinates 
(columns of the generator matrix), followed by multiplication of coordinates by 
nonzero elements from GF(4), followed by possible conjugation of coordinates. 
For a code of length n, there is a total of 6™n! such maps. The 6 possible 
transformations given by scaling and conjugation of a coordinate are equivalent 
to the 6 permutations of the elements {1,lo,co 2 } in the coordinate. A map that 
maps C to C is called an automorphism of C. All automorphisms of C make up an 
automorphism group, denoted Aut(C). The number of distinct codes equivalent 
to C is then given by | Autre) I • ^ v summing the sizes of all equivalence classes, 
we find the total number of distinct codes of length n, denoted T n . It was shown 
by H5hn that T n is also given by the mass formula, 



T„=n(2 i +i)=£' 



6 n n\ 



Aut(Cj 



(5) 



where the sum is over all equivalence classes. Similarly, the total number of 
distinct Type II codes of length n is given by 



n-l 



6™n! 



1 1 Author 



(0) 



where the sum is over the equivalence classes of Type II codes. By assuming 
that | Aut(Cj)| = 1 for all j in Eq. (0, we get a lower bound on t n , the number 
of inequivalent codes of length n. 



tn > 



niu(2' + i) 

6™n! 



(7) 



A similar bound on t 1 ^ can be derived from Eq. JJjJ. 

We can use the computational algebra system Magma to find the auto- 
morphism group of a code. Since, at this time, Magma has no explicit function 
for calculating the automorphism group of an additive code, we use the following 
method, described by Calderbank et al. We map the (n, 2 k ) additive code C 
over GF(4) to the [3n, k) binary linear code /3(C) by applying the map i— > 000, 
1 i— * Oil, co i — > 101, to 2 i — * 110 to each generator of C. We then use Magma to 
find Aut(/3(C)) n Aut(/3(GF(4)")), which will be isomorphic to Aut(C). 

If we are given t n inequivalent codes of length n, i.e., one code from each 
equivalence class, it is relatively easy to calculate the automorphism group size 
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of each code, as described above, and verify that the mass formula defined by 
Eq. JSJ gives the correct value. But to actually find a set of t n inequivalent codes, 
or just the value of t n , is a hard problem. All self-dual additive codes over GF(4) 
of length n were first classified, up to equivalence, by Calderbank et al. Q for 
n < 5 and by H5hn for n < 7. Hohn also classified all Type II codes of 
length 8. Using a different terminolo gy, the codes of length n were implicitly 
classified by Hein, Eisert, and Briegel [14j for n < 7 and by Glynn et al. [12( for 
n < 9. These classifications were not verified using the mass formula defined 
by Eq. Q. Gaborit et al. [1GI3 

have classified all extremal codes of length 8, 
9, and 11, and all extremal Type II codes of length 12. Bachoc and Gaborit 
classified all extremal Type II codes of length 10, and they also showed that 
there are at least 490 extremal Type II codes of length 14 and gave a partial 
result on the unicity of the extremal Type II code of length 18. A review of the 
current status of the classification of various types of self-dual codes is given by 
Huffman 0. 

In this paper, we will give a complete classification of all codes of length up 
to 12, and all extremal Type II codes of length 14. But first, in Section [21 we 
introduce isotropic systems and show that they correspond to self-dual additive 
codes over GF(4). It is known that isotropic systems can be represented by 
graphs. In Section we define graph codes. Theorem El shows that every code 
can be represented by a graph. This gives us a much smaller set of objects to 
work with. In Section 0] we introduce local complementation, and Theorem 1 121 
states that two codes are equivalent if and only if the corresponding graphs 
are related via local complementations and graph isomorphism. This implies 
that classifying codes up to equivalence is essentially the same as classifying 
orbits of graphs under local complementation. We describe an algorithm for 
generating such graph orbits in Section This algorithm was used to classify 
all codes of length up to 12. We show that Type II codes correspond to a special 
class of graphs and use this fact to classify all extremal Type II codes of length 
14. Finally, we determine that the smallest Type I and Type II codes with 
trivial automorphism group have length 9 and 12, respectively. In Section we 
conclude and mention some other results. 

2 Isotropic Systems 

We define a mapping : GF(4) -> GF(2) 2 by <j>{x) = (TrfW 2 ), Tr(x)), i.e., 
i — ► (0,0), 1 i — ► (1,0), u> i — ► (0,1) and ui 2 i— > (1,1). The reverse mapping 
(j)- 1 : GF(2) 2 -► GF(4) is given by ^ (a, 6) = a + Lub. Let u € GF(2) 2 ™ be writ- 
ten as u = (a\b) = (ai, 02, ... , a n , b\, 62, ■ ■ ■ , b n ). We extend the mapping <j> : 
GF(4) n -> GF(2) 2 ™ by letting <p(v) = (a\b) where <j>(vi) = {a h bi). Likewise, we 
define ^ : GF(2) 2 " -> GF(4) n by -1 (a|&) = a+cob. We define the symplectic 
inner product of (a\b), (a'\b') g GF(2) 2 ™ as ((o|6), (o'|6')> = a ■ b' + b ■ a'. A 
subset X C GF(2) 2 ™ is called totally isotropic if (it, v) = for all it, v G X. 

Definition 2. A totally isotropic linear subspace of GF(2) 2n with dimension n 
defines an isotropic system . An isotropic system can therefore be defined by 
the row space of a full rank nx2n binary matrix (A\B), where AB T + BA T = 0. 

Theorem 3. Every self- dual additive code over GF(4) can be uniquely rep- 
resented as an isotropic system, and every isotropic system can be uniquely 



4 



represented as a self-dual additive code over GF(4). 



Proof. Let C C GF(4) n be a self-dual additive code. Map Ctolc GF(2) 2n by 
mapping each codeword u G C to 4>(u) = (a\b) G GF(2) 2 ™. X must then be a 
linear subspace of GF(2) 2 ™ with dimension n. (a\b), (a'\b') G T are orthogonal 
with respect to the symplectic inner product if and only if <f)~ 1 (a | b) , cj>~ 1 (a' | b' ) G 
C are orthogonal with respect to the Hcrmitian trace inner product, because 

(p- 1 (a\b)*(l)- 1 (a'\b') 



= Tr^-^olft) -^-^o'lb')) 
= Ti({a + Lub) ■ (a' +Ub')) 

= {a ■ a') Tr(l) + (a • b') Tr(w) + (6 • a') Tr(w) + (b ■ b') Tr(l) 
= a ■ b' + b ■ a'. 



Since C is self-dual, u * v = for all u,teC, and 2" must therefore be totally 
isotropic. It follows that X defines an isotropic system. Likewise, the reverse 
mapping from an isotropic system to a subset of GF(4)™ will always give a 
self-dual additive code over GF(4). □ 

Example 4. The row-space of (A\B) defines an isotropic system, while C = 
A + uoB is a generator matrix of the (6, 2 6 , 4) Hexacodc. 



(A\B) 



/100100 
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10 10 
10 1 
10 1 

\0 1 1 



1 1\ 
10 111 
10 1 
10 111 
110 
1111/ 
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3 Graph Representation 

A graph is a pair G = (V, E) where V is a set of vertices, and C V x 1/ is a 
set of edges. A graph with n vertices can be represented by an n x n adjacency 
matrix T, where 7^ = 1 if {i,j} G E, and 7^ = otherwise. We will only 
consider simple undirected graphs whose adjacency matrices are symmetric with 
all diagonal elements being 0. The neighbourhood of v G V, denoted C V, is 
the set of vertices connected to v by an edge. The number of vertices adjacent 
to v, \N V \, is called the degree of v. The induced subgraph of G on W C V 
contains vertices W and all edges from E whose endpoints are both in W. The 
complement of G is found by replacing E with V x V — E, i.e., the edges in E 
are changed to non-edges, and the non-edges to edges. Two graphs G = (V, E) 
and G' = (V, E') are isomorphic if and only if there exists a permutation it 
of V such that {u,v} G E ^==^> {ir(u), tt(v)} G E' . A pa£/i is a sequence of 
vertices, (vi, v 2 , ■ ■ ■ , vi), such that {v\, v 2 }, {v2,v 3 }, . . . , e E. A graph 

is connected if there is a path from any vertex to any other vertex in the graph. 

Definition 5. A graph code is an additive code over GF(4) that has a generator 
matrix of the form C = T + u>I, where I is the identity matrix and T is the 
adjacency matrix of a simple undirected graph. 
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A graph code is always self-dual, since its generator matrix has full rank over 
t 

GF(2) and CC only contains entries from GF(2) whose traces must be zero. 
This construction for self-dual additive codes over GF(4) has also been used by 
Tonchev 0. 

Theorem 6. Every self-dual additive code over GF(4) is equivalent to a graph 
code. 

Proof. (This proof is due to Van den Nest, Dehaene, and De Moor [25I Eij.l 
We recall that the generator matrix of a self-dual additive code over GF(4) 
corresponds to an n x 2n binary matrix (A\B), such that C = A + loB. The 
row-space of (A\B), denoted X, defines an isotropic system. We must prove 
that X is also generated by (r|/), where J is the identity matrix and T is the 
adjacency matrix of a simple undirected graph. 

The rows of (A|_B) can be replaced by any n independent vectors from X. 
This basis change can be accomplished by (A'\B') = M(A\B), where M is 
an n x n invertible binary matrix. If B is invertible, the solution is simple, 
since B~ 1 (A\B) = (T\I). Note that T will always be a symmetric matrix, since 
r/ T + iT T = 0. If the ith diagonal element of T is 1, it can be set to by 
conjugating column i of T + u)I. 

In the case where B has rank k < n, we can perform a basis change to get 



where B\ is a k x n matrix with full rank, and A\ also has size k x n. Since 
the row-space of (A'\B') is totally isotropic, and B' contains an all-zero row, it 
must be true that A2BJ — 0. A2 must have full rank, and the row space of B\ 
must be the orthogonal complement of the row space of A2 . 

We assume that B\ — (-B11I-B12) where B\\ is a k x k invertible matrix. 
We also write A2 — (A21IA22) where A22 has size (n — k) x (n — k). Assume 
that there exists an x G GF(2) ,l_fc such that A22X 1 - = 0. Then the vector 
v = (0, . . . , 0, x) of length n satisfies A2V 1 - = 0. Since the row space of B\ is 
the orthogonal complement of the row space of A2 , we can write v — yB\ for 
some y € GF(2) fc . We see that yB\\ = 0, and since B\\ has full rank, it must 
therefore be true that y — 0. This means that x = 0, which proves that A22 is 
an invertible matrix. 

Interchanging column i of A' and column i of B' corresponds to multipli- 
cation by to 2 followed by conjugation of the ith column of A' + loB 1 . We can 
therefore swap the ith columns of A 1 and B 1 for k < i < n to get (A"\B"). Since 
-Bn and A22 are invertible, B" must also be an invertible matrix. We then find 
B"- 1 {A"\B") = (r|J), and set all diagonal elements of Y to 0. □ 

Example 7. Let C = A+ujB be the generator matrix of the (6, 2 6 , 4) Hexacode 
given in Example ^ By the method described in the proof of Theorem El we 
find C = r + u)I, which generates an equivalent graph code. T is the adjacency 
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(a) (b) 

Fig. 1: Two Graph Representations of the Hexacode 



matrix of the graph shown in Fig. Ilbl 
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Theorem[(|]was first proved by Bouchet [3j in the context of isotropic systems, 
and later by Schlingemann |20| in terms of quantum stabilizer states. Proofs of 
Theorem El have also been given by Schlingemann and Werner |2lJ, by Grassl, 
Klappenecker, and Rotteler 0], by Glynn et al. 0,0], and by Van den Nest 
et al. [Hill. 

Swapping vertex i and vertex j of a graph with adjacency matrix T can be 
accomplished by exchanging column i and column j of T and then exchanging 
row i and row j of T. We call the resulting matrix V . Exactly the same column 
and row operations map T + ujI to T' + luI. These matrices generate equivalent 
codes. It follows that two codes are equivalent if their corresponding graphs are 
isomorphic. 

We have seen that every graph represents a self-dual additive code over 
GF(4), and that every self-dual additive code over GF(4) can be represented by 
a graph. It follows that we can, without loss of generality, restrict our study 
to codes with generator matrices of the form T + tol, where T are adjacency 
matrices of unlabeled simple undirected graphs. 
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(a) The Graph G 



(b) The LC Image G° 



Fig. 2: Example of Local Complementation 

4 Local Complementation 

Definition 8. Given a graph G = (V, E) and a vertex v € V, let N v C V be the 
neighbourhood of v. Local complementation (LC) on v transforms G into G v . 
To obtain G v , we replace the induced subgraph of G on N v by its complement. 
It is easy to verify that (G v ) v = G. 

Example 9. We will perform local complementation on vertex of the graph 
G, shown in Fig. [23 We see that the neighbourhood of is N = {1, 2, 3} and 
that the induced subgraph on the neighbourhood has edges {1,2} and {1,3}. 
The complement of this subgraph contains the single edge {2, 3}. The resulting 
LC image, G°, is seen in Fig. I2bl 

Example 10. Consider the graph shown in Fig.^| whose corresponding graph 
code is the Hexacode. An LC operation on any vertex of this graph produces 
the graph shown in Fig. llbl An LC operation on the vertex in the centre of the 
graph shown in Fig. llbl gives the same graph, up to isomorphism. LC operations 
on any of the other five vertices produces the graph shown in Fig. Hal 

Theorem 11. Let T be the adjacency matrix of the graph G — (V, E), and T v be 
the adjacency matrix of G v , for any v EV. The codes generated by C = T + ujI 
and C' = T" + ujI are equivalent. 

Proof. We show that C can be transformed into C' by using only operations that 
map a code to an equivalent code. Each row and each column of C correspond 
to a vertex in V . Let N v denote the neighbourhood of v. For all i € N v , add 
row v of C to row i of C . Multiply column v of C by uj and then conjugate the 
same column. Finally, conjugate column i of C, for all i S N v . The resulting 
matrix is C' . □ 

Theorem 12. Two self-dual additive codes over GF(4), C and C , with graph 
representations G and G' , are equivalent if and only if there is a finite sequence 
of not necessarily distinct vertices (v\, i>2, • ■ • ,Vi), such that (((G Vl ) V2 )'") Vi is 
isomorphic to G' . 

Sketch of proof . Let T be the adjacency matrix of G, and let Cq be the code 
generated by T + ojI. Likewise, let T' be the adjacency matrix of G' , and let C' G 
be the code generated by T' + uil. If the codewords of C are mapped onto the 
codewords of C by one of the 6™n! combinations of coordinate permutations, 
coordinate scalings, and coordinate conjugations, then there must also be a 
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transformation from this set that maps the codewords of Co onto the codewords 
of C'q. Consequently, we only need to consider those transformations that map 
a graph code to another graph code. The codes obtained by the n\ possible 
permutations of coordinates correspond to graph isomorphisms. 

Let C = r + uil be transformed into C = A + uiB by coordinate scalings 
and conjugations. Then C is a graph code if and only if B is invertible and 
all diagonal elements of B~ 1 A are zero. It is easy to verify that conjugation 
of column i of C has no effect on B, but flips the value of the ith diagonal 
element of B~ x A. Given a combination of column scalings on C such that 
the resulting B is invertible, there must therefore be a unique combination of 
column conjugations on C such that the resulting B~ 1 A has zero diagonal. We 
must therefore show that any combination of column scalings on C that give an 
invertible B can be performed by a sequence of LC operations on G. 

Multiplying column i of C by uj 2 replaces column i of I with column i of 
r. Multiplying column i of C by uj adds column i of T to column i of I. It is 
then possible to show which of the 3™ possible scalings do not give an invertible 
B. A vertex v of G corresponds to a column of Y. The neighbourhood of v, 
N v , corresponds to a set of columns of T. We know from Theorem II II that an 
LC operation on vertex i of G corresponds to a scaling of column i of C by 
lu followed by conjugation of column i and all columns in iV,-. Observe that 
conjugating a coordinate followed by a scaling by uj is equivalent to scaling by 
uj 2 followed by conjugation. Note in particular that the local complementations 
((G i yy, where i and j are adjacent vertices, are equivalent to scaling both 
column i and column j of C by uj 2 . It can be verified that any combination of 
column scalings that map a graph code to a graph code can be implemented 
as a sequence of LC operations. The exact algorithm for finding this sequence 
of LC operations is quite involved, and we refer to the proof by Van den Nest 
et al. [Hill for the details. □ 

Bouchet || first proved Theorem ^] in terms of isotropic systems. The 
same result was discovered by Van den Nest et al. [ijj [2(| in terms of quantum 
stabilizer states, and by Glynn et al. [lllll2| using finite geometry. 



5 Classification 

Definition 13. The LC orbit of a graph G is the set of all unlabeled graphs 
that can be obtained by performing any sequence of LC operations on G. 

It follows from Theorem El that two self-dual additive codes over GF(4) are 
equivalent if and only if their graph representations are in the same LC orbit. 
As an example, the two graphs shown in Fig. II al and Fig. llbl make up a complete 
LC orbit, and are thus the only possible graph representations of the Hexacode. 
The LC orbit of a graph can easily be generated by a recursive algorithm. We 
have used the program nauty [18| to check for graph isomorphism. 

Let G n be the set of all unlabeled simple undirected connected graphs on 
n vertices. Connected graphs correspond to indecomposable codes. A code 
is decomposable if it can be written as the direct sum of two smaller codes. 
For example, let C be an (n, 2", d) code and C an (n',2 n , d') code. The di- 
rect sum, C © C = {u\\v | u £ C,v G C'}, where || means concatenation, is 
an (n + n',2 n+n , min{d, d'}) code. It follows that all decomposable codes of 
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length n can be classified easily once all indecomposable codes of length less 
than n are known. 

The set of all distinct LC orbits of connected graphs on n vertices is a 
partitioning of G n into i n disjoint sets. i n is also the number of indecomposable 
self-dual additive codes over GF(4) of length n, up to equivalence. Let L n 
be a set containing one representative from each LC orbit of connected graphs 
on n vertices. We have devised several algorithms Q for finding such sets of 
representatives. The simplest approach is to start with the set G n and generate 
LC orbits of its members until we have a partitioning of G n . The following more 
efficient technique was described by Glynn et al. [12| . Let the 2™ — 1 extensions 
of a graph on n vertices be formed by adding a new vertex and joining it to all 
possible combinations of at least one of the old vertices. The set E n , containing 
i ra _i(2 n_1 — 1) graphs, is formed by making all possible extensions of all graphs 
in L n -i- 

Theorem 14. L n C E n , i.e., the set E n will contain at least one representative 
from each LC orbit of connected graphs on n vertices. 

Proof. Let G — (V, E) 6 G n , and choose any subset W C V of n— 1 vertices. By 
doing LC operations on vertices in W, we can transform the induced subgraph 
of G on W into one of the graphs in L n -\ that were extended when E n was 
constructed. It follows that for all G £ G n , some graph in the LC orbit of G 
must be part of E n . □ 

The set E n will be much smaller than G n , so it will be more efficient to 
search for a set of LC orbit representatives within E n . It is also desirable to 
partition the set E n such that graphs from two different partitions are guar- 
anteed to belong to different LC orbits. We can then consider each partition 
independently, which reduces the amount of memory required and allow for par- 
allel processing. To do this, we must have some property that is invariant over 
the LC orbit and that can be calculated quickly. 

The special form of the generator matrix of a graph code makes it easier to 
find the number of codewords of weight i < n. If C is generated by C = T + uil, 
then any codeword formed by adding i rows of C must have weight at least i. 
This means that we can find the partial weight distribution of C, (Aq, Ax, . . . , Aj), 
for some j < n, by only considering codewords formed by adding j or fewer 
rows of C. We calculate the partial weight distribution, for a suitable choice 
of j, of all codes corresponding to graphs in E n . Codes with different partial 
weight distribution can never be equivalent, so we partition E n such that graphs 
corresponding to codes with the same partial weight distribution are always in 
the same partition. 

Using the described techniques, and a parallel cluster computer, we were 
able to classify all self-dual additive codes over GF(4) of length up to 12. The 
results have been verified by checking that the sizes of all LC orbits add up to 
the number of graphs in G n . The sizes of the automorphism groups of all codes 
have also been calculated, and it has been verified that that the mass formulas 
defined by Eq. (JSJ and Eq. © give the correct values. Table Hgi ves the values 
of i ni the number of distinct LC orbits of connected graphs on n vertices, which 
is also the number of inequivalcnt indecomposable codes of length n. The table 
also gives the values of i 1 ^, the number of indecomposable Type II codes. The 
total number of inequivalent codes of length n, t n , and the total number of 
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Table 1: Number of Indecomposable (i„) and Indecomposable Type II (i 1 ^) Codes of 
Length n 



nl23456 7 8 9 10 11 12 



i n 1 1 1 2 4 11 26 101 440 3,132 40,457 1,274,068 
i* 1 1 1 4 14 103 2,926 



Table 2: Total Number (t„) and Number of Type II (t„) Codes of Length n 



nl2345 6 7 8 9 10 11 12 



t n 1 2 3 6 11 26 59 182 675 3,990 45,144 1,323,363 
t^ 1 1 2 6 21 128 3,079 



Type II codes of length n, t 1 ^, are shown in Table The numbers t n are easily 
derived from the numbers i n by using the Euler transform j23j | , 



d\ 

h = ci 



l( ^ 

V k=l 



The numbers ijj 1 are similarly derived from i 1 ^. The values of i n and t n can 
be found as sequences A090899 and A094927 in The On-Line Encyclopedia of 
Integer Sequences [12] • Table 01 and Table 0] list by minimum distance the 
numbers of indecomposable codes and the total numbers of codes. 1 Tableland 
Table similarly list the numbers of Type II codes by minimum distance. The 
numbers of Type I codes can be obtained by subtracting the numbers of Type II 
codes from the total numbers. The number of distinct weight enumerators of all 
codes of length n and minimum distance d can be found in Table There are 
obviously too many codes to give a complete list here, but a database containing 
one representative from each equivalence class, with information about weight 
enumerators, automorphism groups, etc., is available on-line [6j. 

Our results give a complete classification of the extremal Type I (10, 2 10 , 4) 
and (12, 2 12 , 5) codes. These classifications were previously unknown. The 101 
extremal Type I (10, 2 10 , 4) codes have 6 distinct weight enumerators, 

Wio,l(x,y) = x 10 + 10x 6 y 4 + 72x 5 y 5 + 160x 4 y e + 2A0x 3 y 7 + 285x 2 y 8 + 200xy 9 + 56y 10 , 

W 10i2 (x,y) = x 10 + lAx 6 y 4 + 64a; 5 y 5 + 156x 4 y li + 256a; 3 y 7 + 281x 2 y 8 + 192xy 9 +6Ch/ 10 , 

Wl ,3 (x, y) = x 10 + 18x 6 y 4 + 56a; 5 j/ 5 + 152a; 4 y 6 + 272x 3 y 7 + 277x 2 y 8 + 184xy 9 + 64y 10 , 

W WA (x,y) = x 10 + 22x 6 y 4 + 48a; 5 j/ 5 + 148a; 4 y 6 + 288a; 3 y 7 + 27Zx 2 y % + 176a^ 9 + 68y 10 , 

Wx 0)5 (x,y) = x 10 + 26x 6 y 4 + 40a; 5 j/ 5 + 144a; 4 y 6 + 304a; 3 y 7 + 269x 2 y s + W8xy 9 + 72y 10 , 

Wio,e(x,y) = x 10 + 30x 6 y 4 + 32a; 5 j/ 5 + 140a; 4 y 6 + 320a; 3 y 7 + 265x 2 y a + WOxy 9 + 76y 10 . 

Table|Hllists the number of such codes by weight enumerator and automorphism 
group size. The 63 extremal Type I (12, 2 12 , 5) codes have 2 distinct weight 

^^Note that some authors ll(t Il5l give 3 as the total number of self-dual (7, 2 7 ,3) codes. 
The correct number is 4. 
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Table 3: Number of Indecomposable Codes of Length n and Minimum Distance d 



d\n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


2 


1 


1 


2 


3 


9 


22 


85 


363 


2,436 


26,750 


611,036 


3 








1 


1 


4 


11 


69 


576 


11,200 


467,513 


4 










1 




5 


8 


120 


2,506 


195,455 


5 




















1 


63 


6 






















1 


All 


1 


1 


2 


4 


11 


26 


101 


440 


3,132 


40,457 


1,274,068 





Table 4: 


Tote 


il 


Numbc 


sr of C 


odes 


; of Length n 


and Minimum Dists 


mce d 


d\n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


1 


1 


1 


2 


3 


6 


11 


26 


59 


182 


675 


3,990 


45,144 


2 




1 


1 


3 


4 


13 


29 


107 


416 


2,618 


27,445 


615,180 


3 










1 


1 


4 


11 


69 


577 


11,202 


467,519 


4 












1 




5 


8 


120 


2,506 


195,456 


5 






















1 


63 


6 
























1 


All 


1 


2 


3 


6 


11 


26 


59 


182 


675 


3,990 


45,144 


1,323,363 



Table 5: Number of Indecomposable Type II Codes of Length n and Minimum Dis- 
tance d 



d\n 


2 


4 


6 


8 


10 


12 


14 


2 


1 


1 


3 


11 


84 


2,133 


? 


4 






1 


3 


19 


792 


? 


6 












1 


1,020 


Total 


1 


1 


4 


14 


103 


2,926 


? 



Table 6: Total Number of Type II Codes of Length n and Minimum Distance d 



d\n 


2 


4 


6 


8 


10 


12 


14 


2 


1 


2 


5 


18 


109 


2,285 


? 


4 






1 


3 


19 


793 


? 


6 












1 


1,020 


Total 


1 


2 


6 


21 


128 


3,079 


> 1,727,942 
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Table 7: Number of Distinct Weight Enumerators of All Codes of Length n and 
Minimum Distance d 



d\n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


1 


1 


1 


2 


3 


5 


10 


23 


46 


116 


320 


909 


3,312 


2 




1 


1 


2 


4 


11 


21 


64 


187 


549 


2,249 


11,419 


3 










1 


1 


2 


4 


15 


33 


125 


625 


4 












1 




2 


2 


7 


28 


178 


5 






















1 


2 


6 
























1 


All 


1 


2 


3 


5 


10 


23 


46 


116 


320 


909 


3,312 


15,537 



enumerators, 

Wi2,i{x, y) = x 12 + 40a; V + 212a; 6 y 6 + 424a; V + 725x 4 y 8 + 1080a; V+ 

980a;V + 504a;?/ 1 + 130y 12 , 
W 12 , 2 {x, y) = x 12 + 48a; V + 188a; V + 432a; V + 765a; V + 1040a; V+ 

972a;V° + 528a;?/ 1 + 122y 12 . 

TableEllists the number of such codes by weight enumerator and automorphism 
group size. 

By observing that graphs corresponding to Type II codes have a special 
property, we are able to extend our classification to all the 1,020 extremal Type II 
(14, 2 14 , 6) codes. It was previously shown by Bachoc and Gaborit [l| that there 
are at least 490 such codes. 

Theorem 15. Let T be the adjacency matrix of the graph G. The code C 
generated by C = T + ujI is of Type II if and only if G is anti-Eulerian, i.e., if 
all its vertices have odd degree. 

Proof. If C is of Type II, then every row of C must have even weight. It follows 
that every row of T must have odd weight, and therefore correspond to an anti- 
Eulerian graph. Conversely, if all rows of C have even weight, C must be of 
Type II, since the codeword formed by adding any subset of these rows must 
also have even weight. This follows from the fact that for any two codewords 
of a self-dual code, there must be an even number of coordinates where the 
codewords have different non-zero values. □ 

An anti-Eulerian graph is the complement of an Eulerian graph, i.e., a graph 
where all vertices have even degree. It is easy to show that all anti-Eulerian 
graphs must have an even number of vertices, and it follows that all Type II 
codes must have even length. To classify Type II codes of length 14, we proceed 
as follows. We take the set L\ 2 containing 1,274,068 LC orbit representatives 
of graphs on 12 vertices. All these graphs are then extended, but in a slightly 
different way than earlier. To each graph we add one vertex and join it to 
all possible combinations of at least one of the old vertices. To each obtained 
graph we then add a second vertex and join it to those of the 13 other vertices 
that have even degree. (If the result is not a connected anti-Eulerian graph, it 
is rejected.) By an argument similar to Theorem 1141 it can be shown that all 



13 



Table 8: Number of Extremal Type I (10, 2 10 , 4) Codes with Weight Enumerator w 
and Automorphism Group of Size a 



a\w 


Wio,i 


W W ,2 


Wi ,3 


Wi ,4 


W W ,B 


Wio,6 


All 


I 
















2 


2 


9 


7 


2 






20 


4 


5 


9 


7 


1 






22 


6 








1 






2 


8 


I 


4 


3 




1 




9 


12 




1 










1 


16 




1 


6 


5 


3 




16 


32 




2 


2 


1 


2 




9 


40 














1 


48 








3 






4 


64 






2 








2 


128 




2 










2 


192 




1 


2 




1 




4 


256 












2 


2 


320 


1 










1 


2 


384 












1 


1 


3840 












1 


1 


All 


15 


32 


29 


13 


7 


5 


101 



Table 9: Number of Extremal Type I (12, 2 12 , 5) Codes with Weight Enumerator w 
and Automorphism Group of Size a 



a\w 


Wi 2 ,l 




All 


1 




25 


25 


2 




23 


23 


3 




1 


1 


4 


3 


4 


7 


6 


1 


3 


4 


8 




2 


2 


24 




1 


1 


All 


4 


59 


63 
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Table 10: Number of (14, 2 14 , 6) Codes with Automorphism Group of Size a 



a 


i 

i 


02o 


o 

A 


258 


>_> 

6 


2 / 


4 


Q Q 
00 


c 



2 / 


o 



lo 


12 


( 


18 


1 


21 


1 


24 


16 


28 


1 


3G 


1 


48 


1 


84 


1 


168 


1 


2184 


1 


6552 


1 


All 


1020 



graphs corresponding to Type II codes of length 14 must be part of this extended 
set. Classifying all Type II codes of length 14 turned out to be infeasible with 
our computational resources. Even when using partitioning by partial weight 
distribution, the largest partitions were too large to be processed. However, we 
were able to generate the LC orbits of all graphs corresponding to (14, 2 14 , 6) 
codes. Extremal Type II codes have a unique weight enumerator, and the weight 
enumerator of a (14, 2 14 , 6) code must be 

W u (x, y) = x 14 + 273x 8 y e + 2457a; V + 7098a; 4 ?/ 10 + 6006a; V 2 + 549?/ 4 . 

Table ITU1 lists the number of codes by automorphism group size. Note that codes 
with 21, 168, and 2184 automorphisms were previously unknown. Generator 
matrices of the codes are available on-line |6J. 

As mentioned before, the set of self-dual linear codes over GF(4) is a subset 
of the self-dual additive codes of Type II. Note that conjugation of single coordi- 
nates does not preserve the linearity of a code. It was shown by Van den Nest j2{| 
that the code C generated by a matrix of the form T + loI can not be linear. 
However, if there is a linear code equivalent to C, it can be found by conju- 
gating some coordinates. Conjugating coordinates of C is equivalent to setting 
some diagonal elements of T to 1. Let A be a binary diagonal matrix such that 
r + A + luI generates a linear code. Van den Nest [2j| proved that C is equiv- 
alent to a linear code if and only if there exists such a matrix A that satisfies 
T 2 + AT + TA + T + I = 0. A similar result was found by Glynn et al. [l^ . 
Using this method, it is easy to check whether the LC orbit of a given graph 
corresponds to a linear code. However, self-dual linear codes over GF(4) have 
already been classified up to length 16, and we have not found a way to extend 
this result using the graph approach. 
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We remark that if C is a self-dual additive code over GF(4) with generator 
matrix T + ujI, it can be shown that the additive code over Z4 generated by 
2T + I has the same weight distribution as C. It has also been shown |l9j 
that self-dual additive codes over GF(4) can be mapped to isodual binary linear 
codes, i.e., codes that are equivalent to their duals, by the mapping 1— > 00, 
1 1 ► 11, uj 1 ► 01 and uj 2 1— > 10. A code over Z4 and a binary code obtained 
from the same self-dual additive code over GF(4) by these two methods are 
related by the well-known Gray map. There are also severals mappings from 
self-dual additive codes over GF(4) to self-dual and self-orthogonal binary linear 
codes [ISIIilOji 

An interesting problem, posed by Hohn [l^|, is to find the smallest code with 
trivial automorphism group, i.e., automorphism group of size 1. We find that 
there is no such code of length up to 8, but there is a single code of length 9 
with trivial automorphism group. This code has generator matrix 

/w 0000001l\ 
w 1 1 
uj 1 1 

Ioj 1 1 
0100^1001 . 
0OOOU1O1 
001001^01 

1 1 1 OwO 

\ 1 1 1 1 1 uj J 

The smallest Type II codes with trivial automorphism groups have length 12. 
One such code is generated by 

/w 0000001001l\ 
OwOOOOlOlOOl 
OOuOllOOOOOl 
OOOwlOlOOOlO 
OOllwOOlOOOO 
OOlOOwOOllll 

1 1 OwO 1 1 1 ■ 

1 1 w 1 1 1 
01000101^101 
OOOOOllllwlO 
HOlOllOOlwO 

\i 1 1 n 1 1 ow/ 

Table ^2 lists the numbers of Type I and Type II codes with trivial automor- 
phism group by length and minimum distance. Note that for length 12, almost 
half the codes have trivial automorphism group. For high lengths, one can ex- 
pect almost all codes to have trivial automorphism group flfij ] . This implies that 
the bound on t n given by Eq. 10 is tighter for higher n. Observe that in Ta- 
ble ^2 no code of minimum distance less than 3 is listed. It is easy to show that 
all codes with minimum distance 1 or 2 must have nontrivial automorphisms. 
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Table 11: Number of Type I (Type II) Codes of Length n and Minimum Distance d 
with Trivial Automorphism Group 



d\n < 8 


9 


10 


11 


12 




14 


3 


1(0) 


113 (0) 


6,247 (0) 


392,649 (0) 


? 


(0) 


4 




3(0) 


1,180 (0) 


163,982 (102) 


? 


(?) 


5 








25 (0) 


? 


(0) 


6 










? 


(625) 


All (0) 


1(0) 


116 (0) 


7,427 (0) 


556,656 (102) 


? 


(?) 



6 Conclusions 

By using graph representation and equivalence via local complementation, we 
have classified all additive codes over GF(4) of length up to 12 that are self- 
dual with respect to the Hermitian trace inner product. It follows from the 
bound given by Eq. (0) that there are at least 72,573,549 codes of length 13. 
It is not feasible to classify all codes of length 13 using our method and the 
computational resources available to us. We were however able to classify the 
1,020 extremal Type II (14, 2 14 , 6) codes. This was done by exploiting the fact 
that Type II codes correspond to anti-Eulerian graphs. Finally, we showed that 
the smallest Type I and Type II codes with trivial automorphism group have 
length 9 and 12, respectively. 

The graph representation of a self-dual additive code over GF(4) can also give 
information about the properties of the code. Tonchev |24j | showed that strongly 
regular graphs give rise to interesting codes. In particular, codes represented 
by the strongly regular Paley graphs are well-known quadratic residue codes. 
We have shown that many extremal an d op timal codes can be represented by 
nested regular graphs 00]. Glynn et al. |l2j | showed that the minimum distance 
of a code is equal to one plus the minimum vertex degree over all graphs in 
the corresponding LC orbit. We have shown that the LC orbit corresponding 
to a code with high minimum distance only contains graphs with both small 
independent sets and small cliques 0, @ ■ 
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